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C/2 . ABSTRACT 

We present the results of an eight-year long monitoring of the radio emission 
from the Luminous Infrared Galaxy (LIRG) NGC 7469, using 8.4 GHz Very Large 
Array (VLA) observations at 0.3" resolution. Our monitoring shows that the late 
time evolution of the radio supernova SN 2000ft follows a decline very similar to that 
displayed at earlier times of its optically thin phase. The late time radio emission of 
SN 2000ft is therefore still being powered by its interaction with the presupernova 
stellar wind, and not with the interstellar medium (ISM). Indeed, the ram pressure 
f*"*- ■ of the presupernova wind is p w v\, ~ 7.6 x 10 _9 dyncm -2 , at a supernova age of 

' t ss 2127 days, which is significantly larger than the expected pressure of the ISM 

ON . around SN 2000ft. At this age, the SN shock has reached a distance r s h ~ 0.06 pc, 

and our observations are probing the interaction of the SN with dense material that 
was ejected by the presupernova star about 5820 years prior to its explosion. From our 
VLA monitoring, we estimate that the swept-up mass by the supernova shock after 
about six years of expansion is M sw « 0.29 Mq, assuming an average expansion speed 



, of the supernova of 10 km s 1 



We also searched for recently exploded core-collapse supernovae in our VLA im- 
ages. Apart from SN 2000ft (S v sw 1760/iJy at its peak, corresponding to 1.1 x 
10 28 erg s _1 Hz -1 ), we found no evidence for any other radio supernova (RSN) more 
luminous than w 6.0 x 10 26 erg s _1 Hz -1 , which suggests that no other Type Iln SN 
has exploded since 2000 in the circumnuclear starburst of NGC 7469. 

Key words: galaxies: Galaxies: Seyfert: individual: NGC 7469 Starbursts: Super- 
novae: individual SN 2000ft: radio continuum: stars radio continuum: galaxies 



1 INTRODUCTION 

The central kiloparsec region of many nearby Luminous 
Infra- Red Galaxies (LIRGs) shows a distribution of their ra- 
dio emission consisting of a compact 150 pc), high surface 
brightness, central radio source immersed in a low surface 
brightness circumnuclear halo (Condon et al. 1991). While 
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the compact, centrally located radio emission in these galax- 
ies might be generated by a point-like source (AGN), or by 
the combined effect of multiple RSNe and supernova rem- 
nants (SNRs), e.g., the evolved nuclear starburst in Arp 220 
(Parra et al. 2007 and references therein) and M 82 (e.g. 
Muxlow et al. 1994), it seems now well established that in the 
circumnuclear regions of those objects there is an ongoing 
burst of star formation producing core-collapse supernovae 
(CCSNe, Type Ib/c and Type II SNe) at a high rate (e.g. 
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NGC 7469, Colina et al. 2001a,b; Arp 299, Neff et al. 2004; 
NGC 6240, Gallimore & Beswick 2004; Arp 220, Parra et al. 
2007; IRAS 18293-3413, Mattila et al. 2007, Perez- Torres 
et al. 2007; IRAS 17138-1017, Kankare et al. 2008a, Perez- 
Torres et al. 2008, Kankare et al. 2008b). Given a reasonable 
assumption for the initial mass function (IMF) (e.g., Colina 
& Perez-Olea 1992; Smith et al. 1998), the detection of CC- 
SNe can be used to characterize the most important param- 
eter in starbursts, namely the star formation rate (SFR). 
However, the estimation of the CCSN rate in LIRGs is a 
challenging task, as the optical emission of supernovae is 
hampered by background brightness and by large amounts 
of dust in the nuclear starburst environment, and therefore 
those SNe remain undiscovered by optical searches. 

Fortunately, it is possible to directly probe the recent 
star forming activity in LIRGs by means of high-resolution, 
high-sensitivity radio searches for CCSNe, since these wave- 
lengths do not suffer from dust obscuration, and can there- 
fore set strong constraints on the properties of star forma- 
tion in the dust-enshrouded environments encountered in 
LIRGs. Significant radio emission from CCSNe is expected, 
as the interaction of the SN ejecta with the circumstellar 
medium (CSM) gives rise to a high-energy density shell, 
which is Rayleigh- Taylor unstable and drives turbulent mo- 
tions that amplify the existing magnetic field in the pre- 
supernova wind, and efficiently accelerate relativistic elec- 
trons, thus enhancing the emission of synchrotron radiation 
at radio wavelengths (Chevalier 1982). Therefore, starburst 
activity in the circumnuclear regions of LIRGs ensures both 
the presence of a high number of massive stars and a dense 
surrounding medium, so bright radio SNe are expected to 
occur (Chevalier 1982, Chugai 1997). In fact, radio and op- 
tical/IR observations have shown that highly-extinguished 
SNe do exist in the circumnuclear (r < 1 kpc) region of local 
LIRGs (e.g., SN 2000ft in NGC 7469, Colina et al. 2001a,b; 
SN 2004ip in IRAS 18293-3413, Mattila et al. 2007, Perez- 
Torres et al. 2007; SN 2008cs in IRAS 17138-1017, Kankare 
et al. 2008a, Perez- Torres et al. 2008, Kankare et al. 2008b). 

NGC 7469 is a well known barred spiral galaxy lo- 
cated at a distance of 70 Mpc (Sanders et al. 2003; Ho= 
70 km s" 1 Mpc" 1 ; at the distance of NGC 7469, 1 arcsec 
corresponds to a linear size of 333 pc) containing a lumi- 
nous, QSO-like, Seyfert 1 nucleus surrounded by a dusty 
circumnuclear starburst of 5" (1.6 kpc) in diameter (Wilson 
et al. 1991, Malkan et al. 1998, Genzel et al. 1995, Scoville 
et al. 2000, Diaz-Santos et al. 2007). NGC 7469 is also a 
LIRG with Li R (8-1000^m) = 4.5x10" L (obtained from 
Sanders et al. 2003 for our adopted distance to NGC 7469). 
During the course of a Very Large Array (VLA) monitoring 
program aimed at the detection of RSNe in NGC 7469, Col- 
ina et al. (2001a) detected a strong compact radio source in 
the circumnuclear starburst of this galaxy, which was later 
identified with a RSN (Colina et al. 2001b). The RSN was as- 
signed the name SN 2000ft (Colina et al. 2002), and was the 
first RSN discovered in the circumnuclear ring of a Seyfert 1, 
Luminous Infrared Galaxy, and its detection confirmed that 
radio observations are a powerful tool for detecting RSNe up 
to large distances in the local Universe. SN 2000ft is located 
at a projected distance of 600 pc from the QSO-like nucleus 
of NGC 7469, and showed an 8.4 GHz emission peak of 1.1 
xlO 28 erg s" 1 Hz" 1 . The case of SN 2000ft is outstanding, 
since its high brightness has allowed a multi-frequency, long 



term monitoring of its radio emission since its discovery. Al- 
berdi et al. (2006) found, based on the flux density evolution 
of SN 2000ft for the first < 1000 days after its explosion, 
that SN 2000ft shares many of the same properties that are 
common to RSNe optically identified as Type II SNe, despite 
having exploded in the dusty and very dense environment 
of a LIRG circumnuclear region. Finally, we also note that 
the optical glow of SN 2000ft has recently been identified in 
archival HST images taken on 13 May 2000 (Colina et al. 
2007), consistent with the predicted date of the explosion 
from the analysis of its radio light curves. 



2 VLA OBSERVATIONS AND DATA 
ANALYSIS 

We have monitored the radio emission from the LIRG 
NGC 7469 since 1998, using the VLAQ in A configuration at 
8.4 GHz (FWHM~0.23-0.31 arcsec). We observed NGC 7469 
on 8 April 1998, 8 September 1999, 27 October 2000, 8 
February 2002, and 19 June 2003 (see Colina et al. 2001b 
and Alberdi et al. 2006), and on 11 November 2004 and 5 
February 2006 (this paper) , using an effective bandwidth of 
100 MHz and both circular polarizations. Each observing 
run lasted about 4 hr (except those on 8 April 1998 and 
8 September 1999, which lasted about 50 and 40 minutes 
each, with effective integration on-target times of ~40 and 
~30 minutes). A typical run consisted of ~16 min scans on 
the nucleus of NGC 7469, interleaved with ~2.5 min scans 
on the phase and amplitude calibrator 2251+158, and each 
time ending with a ~6 min observation of the quasar 3C 48 
to set the absolute flux density scale. This observing scheme 
resulted in effective integration times on NGC 7469 of «3.0 
hr. We later edited, calibrated, and imaged the VLA data 
corresponding to all epochs by following standard data re- 
duction techniques implemented within the NRAO Astro- 
nomical Imaging Processing System (AIPS). 

The radio maps shown in Fig. [T] were obtained after 
a few iterations of imaging and phase-selfcalibration of the 
NGC 7469 data, setting ROBUST=2 in task IMAGR, and 
the achieved synthesized beam ranged between 0.26" and 
0.31", depending on the u-v coverage at a given observ- 
ing epoch. To facilitate comparisons among the observing 
epochs, the maps have been spatially convolved to obtain 
a same circular Gaussian beam of size 0.31" (see Fig. [TJ, 
which corresponds to the lowest resolution image on 8 April 
1998. The off-source rms for each map is of 25, 25, 15, 17, 20, 
14, and 14 /ijy beam" 1 for the 1998, 1999, 2000, 2002, 2003, 
2004, and 2006 images, respectively, which are in general 
close to the theoretically expected image thermal noise val- 
ues. Figure [5] shows the resulting image of combining the ra- 
dio interferometric data for all epochs (1998 through 2006), 
where we have marked for convenience the positions of the 
nucleus (N), the well-known supernova SN 2000ft, and the 
star forming regions Rl, R2, and R3 (Colina et al. 2001b). 



1 The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agree- 
ment by Associated Universities Inc. 
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Figure 1. Total intensity radio images of NGC 7469 obtained at 8.4 GHz with the VLA in A-configuration, which cover a time baseline 
of almost eight years (from 8 April 1998 till 5 February 2006). The maps have been spatially convolved to obtain the same circular 
Gaussian beam of size 0.31", to be able to compare adequately all the images. The contour levels are drawn at (3,3 y/3,9,...) x25fi Jy 
beam -1 , which is the highest off-source rms of all the maps (corresponding to the 1998 and 1999 epochs). The peaks of brightness are 
of 12.61, 12.40, 12.27, 12.53, 13.01, 12.91, and 13.15 mjy beam -1 for the 1998 through 2006 epochs, respectively, and correspond to the 
nucleus of NGC 7469. The cross corresponds to the position of the peak of SN 2000ft in the 2000 image. 
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3 RESULTS 

3.1 Radio emission from NGC 7469 and its SN 
2000ft 

We summarize the results of our eight-year long VLA-A 
monitoring of NGC 7469 and its radio bright SN 2000ft in 
Table [1] and Figure [T] The average measured 8.4 GHz total 
flux density of the galaxy is ~35 mjy. A significant frac- 
tion of this radio emission (~15-17 mjy, depending on the 
epoch), comes from the nuclear region of the galaxy, with an 
angular size of ~0.5" (~167 pc), which hosts an AGN. We 
also note that our MERLIN (Multi-Element Radio Linked 
Interferometer Network) observations at 5 GHz (Alberdi et 
al. 2006), imply a linear size < 50 pc for the nuclear re- 
gion. Further, recent high-resolution images of this region 
obtained with the European VLBI Network (EVN) show 
that the putative AGN consists of at least four compact 
components (Alberdi et al., in preparation). 

Our 8.4 GHz VLA images clearly show the appear- 
ance of a radio supernova, SN 2000ft, whose extreme bright- 
ness and long-lasting behavior confirmed not only its core- 
collapse nature, but also that it was a Type Iln super- 
nova (Alberdi et al. 2006). The progenitors of those SNe 
are thought to be rather massive stars, about (18—30) Mq 
(Weiler et al. 1990, van Dyk et al. 1993), although we should 
note that there is also recent evidence that some of the Type 
Iln progenitors might be much more massive, e.g. Gal- Yam 
et al. (2007), Trundle et al. (2008). We obtained the flux den- 
sity of SN 2000ft for each image in Figure[T] after subtraction 
of the background emission, i.e., the local, secular, non vari- 
able emission of the galaxy at 8.4 GHz. We estimated the 
local background radio emission in the region around SN 
2000ft by four different procedures: (i) We used the AIPS 
task IMFIT, solving for a zero level component (the back- 
ground) for each observing epoch. This procedure yielded 
background values of 123, 96, 88, 100, and 99 microJy for 
the 2000, 2002, 2003, 2004, and 2006 epochs, respectively. 
The above background values were subtracted from the peak 
values estimated for SN 2000ft, and are the values quoted 
in Table [1] (ii) We also determined the background emis- 
sion at the position of the supernova as the difference in the 
flux density of two images, the first one obtained by using a 
purely naturally weighting scheme (UVWTFN='NA'; RO- 
BUST^), and the second one obtained by using a uniform 
weighting scheme (UVWTFN=' '; ROBUST=-5), where in 
addition we dropped the inner 35 kilo-A of the u-v range, 
essentially removing all extended emission from the data. 
The former image yields the lowest off-source r.m.s. in the 
images, and is most sensitive to the extended, secular back- 
ground emission, while the second is essentially free of back- 
ground emission, at the expense of having a higher off-source 
rms. The estimates for the secular background emission ob- 
tained in this way coincided with those obtained by following 
the first procedure within 2 /ijy/b. As a measure of the un- 
certainty in the locally subtracted background emission we 
took the standard deviation of the above values (13/iJy/b). 
We conservatively added this uncertainty linearly, rather 
than in quadrature, to the total uncertainty in the flux den- 
sity measurements for SN 2000ft, which are listed in Table [T] 
and plotted in Figure[3] Those values show that the late time 
radio emission of SN 2000ft is well characterized by a single 
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Figure 2. Total intensity radio image of NGC 7469 obtained at 
8.4 GHz with the VLA in A-configuration, using all u — v data for 
epochs 1998 through 2006, whose individual images are shown in 
Figure [T] The radio image has been obtained by convolving the 
data with a circular Gaussian beam of size 0.31". The contour 
levels are drawn at (3,3 V3,9,...) X the off- source rms of the image, 
which is 10 /tjy beam -1 . This is the deepest 8.4 GHz image of 
NGC 7469 ever produced. We indicate the positions of the nucleus 
(N), the RSN SN 2000ft, and the star forming regions Rl, R2, and 
R3. 

index for the power-law time decay, even six years after its 
explosion (S„ tx t ' -0.07 ■ se e Sect. 13.21 for details). 

We note here that the flux density values reported in 
Table [T] for SN 2000ft are somewhat different from those 
presented in Alberdi et al. (2006). Their values were deter- 
mined from images with different synthesized beams, used 
different weighting imaging schemes, applied amplitude self- 
cal after a few iterations of phase-selfcalibration, and used 
different field and cell sizes for each of the images presented 
in their paper. Here, we followed a homogeneous data re- 
duction scheme. Specifically, we used the same field and cell 
sizes for all images, made a standard imaging scheme, using 
natural weighting and applying phase-selfcalibration only. 
The final run of IMAGR was done using a circular Gaussian 
beam of size 0.31 arcsec for all epochs. Finally, we took into 
account the subtraction of the background radio emission. 
All this explains the differences in the flux densities reported 
in those papers 0, but do not affect, or modify, any of the 
results or conclusions obtained in either paper. 

3.2 Radio light curve of SN 2000ft 

We show in Figure[3]the flux density evolution of SN 2000ft. 
The data plotted correspond to our new flux density esti- 

2 The only exception is the flux density reported for the discovery 
of SN 2000ft, which is mistakenly given as 1.60 mjy in Alberdi et 
al. (2006), while the value that should have been reported is 1.70 
mJy. 
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Table 1. 8.4 GHz positions and flux densities of compact sources in the (circum)nuclear regions of NGC 7469. The uncertainty in the 
flux density determination results from adding in quadrature the off-source rms of each map and a 2% of the local maximum and, linearly, 
the standard deviation of the locally subtracted background emission. 
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Figure 3. Radio light curve of SN 2000ft at 8.4 GHz (filled 
squares), 5.0 GHz (filled circles), and 1.6 GHz (upper limits; open 
triangles), along with the resulting light curve fits (8.4 GHz, solid 
line; 5.0 GHz, dotted line; 1.4 GHz, dashed-dotted line) of the pa- 
rameterized model described in Section l3.2l The explosion date is 
fixed to 10 May 2000, three days before its serendipitous optical 
discovery (Colina et al. 2007). The 8.4 GHz data are from this 
work, while data at other wavelengths were taken from Alberdi 
et al. (2006). Note that the late time data (t > 1000 days) fol- 
lows genuinely the same steep, power-law decline found at earlier 
times. 



mates (corrected for the galaxy background emission), as 
well as our late time (t > 1000 days) radio data at 8.4 GHz 
(see Table[T} , while data at other wavelengths are taken from 
Alberdi et al. (2006). (Note that since the data are other 
frequencies were obtained with MERLIN, which has a reso- 
lution of ~40 mas at 5 GHz, the putative background radio 
emission was resolved out). Following Alberdi et al. (2006), 
we fitted the light curves in terms of the "mini-shell" model 
(Chevalier 1982), including modifications by Weiler et al. 
(1990). In particular, Alberdi et al. (2006) found that both 
internal opacity (consisting of synchrotron self-absorption 
and mixed thermal free-free absorption/nonthermal absorp- 
tion) and clumpy opacity contributions were not required to 
fit the data. Alberdi et al. (2006) did find, however, that a 
distant foreground absorber was needed to account for the 
non-detection of SN 2000ft at 1.6 GHz in any of the observ- 
ing epochs, so we have also fitted for such an absorber. In 
addition, we also use the fact that SN 2000ft has been re- 
cently detected at optical wavelengths on 13 May 2000 (Col- 



ina et al. 2007), based on an analysis of HST archival data. 
We have fixed here the explosion date to 10 May 2000, so 
as to make it compatible with the optical detection. There- 
fore, we fitted our multi-wavelength observations using the 
following expression (e.g. Weiler et al. 2002): 

S v = Ki v% t{ e~ T ^ (1) 

where 5V is the observed flux density, in mjy; ^5 is the ob- 
serving frequency, in units of 5 GHz; t\ is the time since 
the explosion date, in days; and T ex t is the external optical 
depth, assumed to arise from purely thermal, ionized hydro- 
gen that absorb with frequency dependence v^ 2 ' 1 . In turn, 

Text = TCSM + Tdist; TCSM ~ K% V^ 2 ' 1 tf , T diat = KiU~ 2A , 

where K\, K2, and K4 correspond to the flux density (Ki) 
and uniform optical depths {Ki, K4) at 5 GHz, one day 
after the explosion, respectively. The parameter 5 accounts 
for the time dependence of the optical depth for the local 
uniform media. 

Our best fit is yielded by the following set of param- 
eters: spectral index, a = —1.27; power-law time decay, 
= -2.02; Kx = 4.45 x 10 5 mjy; K 2 = 1.67 x 10 7 ; adopted 
power-law time decay of the optical depth, 8 = —2.94; 
foreground uniform optical depth r// (most likely associ- 
ated with an H II region, as found by Alberdi et al. 2006), 
K A > 0.17. 



3.3 A systematic search for radio supernovae in 
NGC 7469 

RSN SN 2000ft was extraordinarily bright, and therefore was 
easy to detect and confirm as such in our images. However, 
RSNe (RSNe) are on average intrinsically fainter than SN 
2000ft, and hence several RSNe could have exploded since 
1998 but gone unnoticed by simple visual inspection of our 
VLA radio images of NGC 7469, especially since there is a 
non-negligible amount (~ lOO^iJy) of secular, non-variable, 
galactic radio emission. We therefore reverted to the use of 
other methods that could hint to the explosion of dimmer 
supernovae during this period. 

In particular, for the purpose of searching for RSNe 
fainter than SN 2000ft, we used all of our seven VLA- A ob- 
serving epochs, including two epochs prior to the discovery 
of SN 2000ft, and re-analyzed them in a similar, homoge- 
neous manner, including editing, calibrating, and imaging 
the data following standard procedures within AIPS. We 
then used a slightly modified version (see Melinder et al. 
2008) of the Optimal Image Subtraction (Alard & Lupton 
1998; Alard 2000) package ISIS 2.2. The method matches 
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the PSF, intensities, and background of the better "seeing" 
image, to the aligned poorer "seeing" images from the dif- 
ferent epochs before the subtraction. The method has been 
used successfully to detect supernovae within the highly ob- 
scured nuclear regions of luminous infrared and starburst 
galaxies in near-IR, e.g. the discovery of SN 2004ip (Mat- 
tila et al. 2007). However, to our knowledge this is the first 
time it is applied to search for supernovae in radio images. 
The radio nucleus of the host galaxy was selected as the im- 
age region to derive the convolution kernel. The epoch from 
2004 was used as the reference frame. No significant effect 
on the image subtraction quality was found when chang- 
ing the kernel basis function width parameters. Flux mea- 
surements were performed on the subtracted images using 
IRAF and GAIA software packages. Due to the incomplete 
(u-v) coverage and the use of the CLEAN deconvolution al- 
gorithm, the noise pattern of the radio images is complex 
and the background is not smooth in the subtracted images 
and contains a large number of false positive and negative 
"sources" . Therefore, to determine our SN detection thresh- 
old we measured the flux in 44 circular apertures around 
the nucleus in the subtracted images using an aperture of 
60 pixels radius and a sky annulus with an inner and outer 
radius of 90 and 120 pixels, respectively (the pixel size was 
of 0.008 arcsec). The standard deviation of these measured 
values was used as a la limit, and corresponded to 20/^Jy/b. 
No additional transient sources besides SN 2000ft above a 
3a detection threshold were apparent in the data set. We 
repeated the above procedure, taking each time a different 
reference epoch, but the results did not change: no addi- 
tional RSNe were detected. 

We also searched for intrinsic variability in regions Rl, 
R2, R3, and N. Since regions Rl, R2, and R3 are very bright 
and show spectral indices typical of star forming regions (Al- 
berdi et al. 2006) , significant variability in their flux density 
well above statistical fluctuations must be due to an intrinsic 
variation caused by, most likely, a supernova. We therefore 
extracted the 8.4 GHz peaks of brightness from the Sy 1 
nucleus of the galaxy (N), and from the star forming re- 
gions Rl, R2, and R3. Those values are also listed in Table 
[1] However, given the uncertainties in the total flux density 
values of those regions, the variations are compatible with 
purely statistical fluctuations. Hence, there is no evidence of 
significant variability that could point towards a supernova. 



4 DISCUSSION 

The main goal of our VLA observations was twofold. First, 
we aimed at monitoring the radio evolution of SN 2000ft, 
especially at late times, and we discuss in detail this as- 
pect in Section 14.21 Second, we searched for recently ex- 
ploded core-collapse supernovae (CCSNe, i.e., Type Ib/c 
and Type II) in the circumnuclear starburst of NGC 7469, 
aimed at establishing its CCSN rate, independently of mod- 
els. Since we obtained in general very high-sensitivity (off- 
source rms < 25/iJy/beam) radio images of NGC 7469 at 
high- resolution (FWHM«0.30"), this could have permitted 
the discovery of new RSNe in the eight years covered by our 
observations. We showed in Section [3] that, apart from SN 
2000ft, there is no evidence for any other bright RSN in the 
(circum) nuclear region of NGC 7469 since 1998. We discuss 



the implications of these results in terms of CCSN rates and 
starburst models in Section \4. 31 

4.1 The mass-loss rate of SN 2000ft 

The new fit to the SN 2000ft light curve (see SectionEOland 
Figure [3]l shows that the best fit parameters are somewhat 
different from those published in Alberdi et al. (2006), due 
to the addition of the new data at 8.4 GHz and also because 
the optical detection of SN 2000ft on 13 May 2000 (Colina et 
al. 2007) implies that the explosion happened before, even if 
close to, the date of the optical detection. While the resulting 
fits do not change any of the facts described in Alberdi et 
al. (2006), one relevant parameter for the interaction of SN 
2000ft with its surroundings deserves particular attention: 
its mass-loss rate. 

In the standard circumsteUar interaction model for 
RSNe, the synchrotron radio emission from the supernova 
is partially free-free absorbed by ionized gas outside the for- 
ward supernova shock front. In this scenario, the supernova 
goes from an optically thick regime to an optically thin one 
when the peak optical depth equals r/^ = f3/5 (this can 
be found by computing the time derivative in Equation [1]) . 
In our case, we have r/^(peak)=0.69. This optical depth is 
attained in the radio regime at an age 

/ ■ \ 2/3 

tf f « 326 I -=± ) T~\ /2 V,- 1 days, (2) 

where M_4 is the mass-loss rate in units of 10~ 4 Mq yr -1 , 
T CS 2 is the circumsteUar temperature in units of 2 x 10 4 K, 
Va is the supernova shock velocity in units of 10 4 km s" 1 , 
v w i is the presupernova wind velocity in units of 10 km s" 1 , 
and f8.45 is the observing frequency in units of 8.45 GHz.. 
Since we are interested in the mass-loss rate, we can write 
the previous equation as follows: 

M « 1.66 x 10" s t% 2 T^V* /2 v wl v s . i5 Mq yr" 1 (3) 

where tioo is the time to reach optical depth unity to free-free 
absorption at a given frequency, in units of 100 days. There- 
fore, we can use our fitted light curves to directly measure 
the time at which SN 2000ft reached rl* ~ 0.69, which cor- 
responds to the emission peak at a given frequency. In par- 
ticular, we find that r i! m 0.69 at «200 days and ~300 days 
at 8.45 and 5.0 GHz, respectively. Substituting those values 
in Equation|3l we find that M < (4.7-5.1) x 10" 5 Mq yr" 1 , 
typical of a red supergiant progenitor star, and where the in- 
equality takes into account that, if synchrotron-self absorp- 
tion plays some role around the observed peak of emission, 
then the obtained mass- loss rate is actually an upper limit. 

4.2 Interaction around SN 2000ft 

It has been proposed that CCSNe exploding in dense en- 
vironments, like those encountered in the circumnuclear re- 
gions of starburst galaxies, can result in very luminous RSNe 
(Chevalier 1982, Chugai 1997). An apparent confirmation of 
such a theory came from the discovery of the radio bright 
SN 2000ft (Colina et al. 2001a,b) in the circumnuclear star- 
burst ring (~600 pc from the nucleus) of NGC 7469. How- 
ever, the VLA radio monitoring of SN 2000ft for the first 
~1100 days after its explosion showed that this RSN shared 
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essentially the same properties that are common to radio 
SNe identified as Type II SNe (Weiler et al. 2002), despite 
having exploded in the dusty and very dense environment of 
the circumnuclear region of NGC 7469. Indeed, Colina et al. 
(2007) have recently identified the optical glow of SN 2000ft 
using archival HST images taken on 13 May 2000, and note 
that the supernova appears to have exploded in a region of 
increased local extinction (~4.2 magnitudes) near the edge 
of a strong lane of dust surrounding the ring (see their color 
map in Figure 1), which confirms that SN 2000ft exploded 
in a dense and dusty surrounding medium. 

Alberdi et al. (2006) also predicted that once the in- 
teraction of the ejecta of SN 2000ft with the interstellar 
medium (ISM) would become relevant, the flux density de- 
cay time scale for SN 2000ft should be longer than for nor- 
mal luminous RSNe, whose radio emission continue to fall 
off very quickly. Such change in the flux density decaying 
rate of the supernova would indicate the termination of the 
presupernova stellar wind, where the pressure of the ISM is 
approximately equal to the ram pressure of the wind, and 
represents the passage of supernova to supernova remnant 
phase. This passage is actually expected to happen earlier 
for supernovae exploding in dense environments-as those 
encountered in starburst galaxies-than in the less dense en- 
vironments of normal galaxies. In fact, the duration of the 
radio SN phase (when the radio emission is governed by in- 
teraction between the ejecta and its circumstellar medium, 
CSM) is limited by the extent of the expanding presuper- 
nova wind, since the latter eventually reaches a radius where 
its ram pressure, p w v^,, equals the external pressure of the 
interstellar medium (ISM), Pism (Chevalier & Fransson 
2001). For a spherically symmetric, steady wind (p w oc r~ 2 ), 
this radius is 

r w n 0.18 M x _!l t&V 1 ^ pc (4) 

where pj = PisAi/k is the ISM pressure in units of 
10 7 cm -3 K, which is the estimated pressure for the central 
region of the starburst in M 82 (Chevalier & Clegg 1985). 
Once it has reached, a core-collapse supernova enters its su- 
pernova remnant (SNR) phase, and its radio emission is no 
longer powered by the interaction with the CSM, but rather 
with the ISM. At this stage, the density encountered by the 
SN would be constant, which would result in a significant 
flattening of the radio emission, i.e., a departure from the 
general, steep decline observed in the early times of the op- 
tically thin phase of SN 2000ft. 



4-2.1 SN 2000ft: Late radio emission still powered by 
interaction with the progenitor's stellar wind 

We followed up the radio evolution of SN 2000ft using the 
VLA at 8.4 GHz, searching for a potential departure of the 
flux density decaying rate, which would signal the passage 
of SN to SNR phase for SN 2000ft. Our observations of the 
late time (t ^ 1000 days) evolution of the 8.4 GHz emission 
of SN 2000ft (see Figure [3]) do not show evidence for any 
departure from the general trend observed for the 8.4 GHz 
emission during its first 1000 days, and implies that the radio 
emission of SN 2000ft is, after more than 2100 days, still 
being powered by the interaction of the supernova with the 
stellar wind of its progenitor star. It then follows that the 



external pressure of the ISM is smaller than the ram pressure 
of the presupernova wind 

PiSM^pu,vl, = M v w /(4:irrw) 

« 5.04 x 10~ 8 M_4 v w j rf 7 2 dyn cm" 2 (5) 

The supernova shock has penetrated into the presuper- 
nova wind up to a distance r s h « 1.84 x 10 17 cm (« 0.06 pc) 
after 2147 days of (assumed) free expansion at a constant 
velocity of 10000 km s" . (We should notice here that the 
radius may be uncertain by as much as a factor of two, 
thus making Pism uncertain by up to a factor of four.) 
Substituting this value, and our estimated mass-loss rate 
of M ~ 4.9 x 10" 5 M yr" 1 for the progenitor star of SN 
2000ft in Equation [5] we then obtain that the ram pres- 
sure of the wind is (p w v^) = 7.6 x 10~ 9 dyn cm -2 , still very 
high to be overcome by Pism- For comparison, Pism ~ 
4 x 10~ 9 dyn cm 2 in the central H II regions of M 82, and 
drops below « 2.1 x 10~ 10 dyncm~ 2 at a distance of 540 pc 
(Chevalier & Clegg 1985). We note here that Parra et al. 
(2007) have reported that some of the most long-lived rem- 
nants in Arp 220 show a rather slow, or even unnoticeable 
flux density decay at 18 cm, as also found by Rovilos et al. 
(2005), and that this implies Pism > 7.6 x 10~ 9 dyncm~ 2 . 
We also note that if the environment where SN 2000ft ex- 
ploded had a similar, or significantly larger, pressure than 
that of the presupernova wind, we should have expected to 
see a significant flattening in the flux density time decay of 
SN 2000ft, which is not observed. 

The number density of the thermal electrons in the stel- 
lar wind at a distance reached by the outer shock can be 
written as 

n w = 3.52 x 10 4 M_4 v~\ r~? cm" 3 (6) 

where rvt = r^/lO 17 cm, and we have assumed solar 
abundances (i.e., the mean molecular weight of the plasma, 
(i = 0.86). Therefore, we obtain n w ~ 5.1 X 10 3 cm -3 at 
an age of t = 2127 days. We stress here that the rele- 
vant physical parameter to determine whether the SN has 
started to enter its SNR phase is not the particle density 
ratio (n m /ni a m), but the ratio of the ram pressure of the 
wind to that of the ISM. Indeed, as we have just shown, the 
thermal particle density in the preshocked wind can be sim- 
ilar to that encountered in the dusty, dense environments of 
LIRGS, yet the radio emission from the supernova be pow- 
ered by its interaction with the CSM. Since we have found 
that Pism /k < 5.5 x 10 7 cm" 3 K, one could have a rather 
dense ISM (msM ~ 10 4 cm" 3 ) at temperatures Tism ~ 10 3 
K, so that this gas would not be able to counterbalance the 
ram pressure of the progenitor wind at the current SN age. 
In this sense, although the conditions of the surrounding 
ISM could be similar to those encountered in the central 
regions of M 82, there is no significant impact on both the 
radio emission and expansion of SN 2000ft, since the ram 
pressure of the presupernova wind is still larger than that 
of the ISM. In order for the ISM to affect noticeably the ra- 
dio emission and expansion of the supernova, and assuming 
the above values of nisM and Tism, SN 2000ft will have to 
expand inside its surrounding ISM up to a radius given by 
Equation|31 which corresponds to r s h — 4.3 x 10 17 cm (= 0.14 
pc), which the supernova will attain at an age of about 13.6 
yr. The SN is therefore expected to follow a similar trend to 
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that observed for its first 5.9 years of expansion until that 
age, at which time it is expected that the SN radio emission 
will start to be dominated by its interaction with the ISM, 
and a significant flattening of the flux should be observed. 

4-2.2 Progenitor mass-loss history and swept-up mass by 
SN 2000ft 

Our late time radio monitoring of SN 2000ft also allowed 
us to shed light on a number of relevant aspects for the 
supernova-CSM interaction, such as the mass-loss history of 
the progenitor star and the swept-up mass by the supernova 
shock, Af sw . If we assume a free expansion for SN 2000ft at 
an average speed of v s = 10 4 km s _1 for the first t ~ 2127 
days (the last epoch of our observations at 8.4 GHz), and 
a standard wind velocity, v w = 10 km s" 1 , it follows that 
the 8.4 GHz emission is currently probing the wind of the 
progenitor star around r = tv s /v w ps 5820 yr prior to its 
explosion. As found in the previous section, the CSM density 
at the corresponding circumstellar wind was n w ~ 5.1 x 
10 3 cm -3 , which is in agreement with expectations for RSG 
supernova progenitors. 

The swept-up mass by the supernova shock is 

P r sh 

M sw =4tt / p w r 2 dr = M rsh ~ r ° M Q (7) 

where we have assumed a standard progenitor wind 
(s — 2; p w oc r~ s ), and r s h, ro, and v m are the current 
shock radius of the supernova (r s h = v s hi), the radius of 
the supernova at shock break-out, and the velocity of the 
stellar wind, respectively. Since r„h S> ro, then Equation [7] 
can be rewritten as 

M sw ~ M r sh /v w = 0.10M.4 v shl v~l fi M (8) 

where v B hi = v s h/W 4 km s -1 , and ti is the time since 
the SN explosion, in years. Since we have M_4 = 0.49, 
v s hi = 1, and £i=5.9, the swept-up mass by the supernova 
shock of SN 2000ft is then M sw ^ 0.29 M Q . 

4.3 The circumnuclear starburst in NGC 7469 

The most relevant parameters of a starburst are its star for- 
mation rate (SFR) and its core-collapse supernova (CCSN) 
rate. The CCSN rate is especially relevant in the study of 
LIRGs, since all these supernovae emit in radio (albeit not 
all are so strong emitters that current radio interferome- 
ters -despite their high sensitivity- are able to detect all of 
them). It is also very important that the (constant) CCSN 
rate, rccsN, can be related to SFR as follows: 

?~ccsn = / <&(m) dm 

Jm SN 

= SFR (Ft ( " >-* ) (9) 

where SFR is the (constant) star formation rate in 
M yr _1 , mi and m„ are the lower and upper mass lim- 
its of the initial mass function (IMF, $ oc m~ a ), and tosn is 
the minimum mass of stars that yield supernovae, assumed 
to be 8 M (e.g., Smartt et al. 2009). Following the formal- 
ism described in Colina & Perez-Olea (1992), the observed 



far-infrared luminosity of a starburst, Lfir, can be used to 
constrain tccsn and SFR. For NGC 7469, which has Lir(8- 
1000/ra) = 4.5 xlO 11 L (Sanders et al. 2003). we obtain a 
constant tccsn > 0.75 SN/yr, and SFR > 40 M yr" 1 , as- 
suming two thirds of the above IR luminosity are emitted by 
its associated circumnuclear starburst (Genzel et al. 1995) , a 
Salpeter IMF (a = 2.35; Salpeter 1955) for m t = 1 M Q and 
m u = 120 M , and that all of the ionizing photons go into 
dust heating. If a fraction around 0.5 of all the photons go 
into dust heating, which is more realistic, then one obtains 
rccsN < 1.5 SN/yr, and SFR < 80 M Q yr -1 . Those values 
are in broad agreement with values previously obtained by 
Wilson et al. (1991), Colina & Perez-Olea (1992), and Gen- 
zel et al. 1995), and also with expectations for core-collapse 
supernova rates based solely in measured far infrared lu- 
minosities, e.g., Mattila & Meikle (2001), which found that 
rccsN » 2.7xl0~ 12 (L F m/L ) yr" 1 , thus implying a CCSN 
rate for the circumnuclear starburst of NGC 7469 of ~0.81 
SN yr"\ for Lir = 3.0 x 10 11 . 

We showed in Section [3] that our 8.4 GHz observations 
indicate the explosion of only one very bright RSN (SN 
2000ft; LP° ak « 1.1 x 10 28 erg s _1 Hz" 1 ) in the circumnu- 
clear starburst of NGC 7469 during the period covered by 
our observations (7.8 years). This implies an upper limit 
for the average rate of very bright RSNe of < 0.13 SN/yr, 
which appears to be far off from the estimated (constant) 
CCSN rate for NGC 7469 of ^(0.81-1.50) SN yr" 1 , based on 
the FIR emission from its circumnuclear starburst. To ex- 
plain this disagreement, it has been usually advocated that 
not all CCSNe result in RSNe. Many CCSNe would then go 
undetected because they will not emit any radio emission. 
However, this is a too simplistic view of CCSNe. Since it is 
the interaction of the ejecta with circumstellar material that 
surrounds CCSN progenitors which gives rise to their radio 
emission, there is no reason to think this physical scenario is 
not valid for all CCSNe, all of which lose significant amounts 
of hydrogen-rich material in the last thousands of years of 
their lives through their stellar winds. Therefore, all CCSNe 
are expected to generate some level of radio emission, al- 
beit not all of them will result in very bright radio emitters, 
which coupled with the limited sensitivity of currently avail- 
able radio interferometers will prevent the detection of faint 
RSNe. 

Indeed, CCSNe display a broad range of radio lumi- 
nosities, spanning several orders of magnitude, with essen- 
tially all detected supernovae showing peak luminosities 
from a few times 10 25 erg s _1 Hz" 1 (all of which would 
go undetected by our observations) up to a few times 
10 28 erg s" 1 Hz" 1 , or more (see e.g., Fig. 5 of Alberdi et al. 
2006). If there was no background emission, RSNe brighter 
than about three times the off-source rms (rms=1.5 x 
10 26 erg s _1 Hz -1 ) could have been directly detected by 
our VLA-A observations. In practice, direct detections of 
RSNe are limited by the background radio emission of the 
galaxy (~ 100/xJy), preventing us from detecting RSNe with 
L R < 6.0xl0 26 erg s" 1 Hz" 1 ). Very bright, long-lived RSNe 
like SN 2000ft are expected to come essentially from Type 
IIL/IIn SNe. Smartt et al. (2009) have found that these rep- 
resent ~6.5% of all CCSNe (assuming a Salpeter IMF), while 
Type IIP/lib SNe -the radio faint ones, with peak luminosi- 
ties of LE eak ~ (5 - 20) x 10 25 erg s" 1 Hz" 1 - are much more 
numerous (~64.1%). (Type Ib/c SNe result often in bright 
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radio events, but they fade away very quickly, typically in 
less than ~100 days, so we do not take them into account, 
as chances of detecting them with our ~1.5 yr monitoring 
of NGC 7469 are meager.) Thus, whatever the CCSN rate, 
about 64.1% of the exploding SNe are expected to become 
relatively radio faint. This would result in a RSN rate of 
« (0.8-1.5) x 0.36=(0. 29-0.54). Since Poisson statistics ap- 
ply to the explosion of SNe, based on a single RSN discovery 
we cannot rule out that a starburst producing massive stars 
at a constant rate is taking place in NGC 7469. 

Another possible explanation within a constant star 
forming scenario for NGC 7469 is that the radio luminos- 
ity function of CCSNe in (U)LIRGs is top-heavy, e.g., if the 
IMF results to be top-heavy, as suggested by, e.g., Klessen 
et al. (2007). In this scenario, we would be witnessing only 
the explosion of very massive stars, which yield also very 
bright radio emitters. We must note, however, that if the 
radio luminosity function of CCSNe is bottom heavy (e.g., 
by mimicking a Salpeter IMF) , then it is very unlikely that 
the only detected event in almost eight years would have 
been SN 2000ft, since similarly bright RSNe should have 
exploded during the 7.8 yrs covered by our VLA observing 
program. 

Finally, we note that while a starburst producing stars 
at a constant rate cannot be excluded, Hubble Space Tele- 
scope multi-wavelength (UV through NIR) imaging and 
if-band ground-based long-slit spectroscopy of NGC 7469 
(Dfaz-Santos et al. 2007) strongly suggest the existence of 
several localized starbursts in NGC 7469, which started at 
different times and in different locations of its circumnu- 
clear ring. According to the observations of Diaz-Santos et 
al. (2007) , the starburst in NGC 7469 consists of many 
massive, young, compact clusters characterized by masses 
between (1-10) xlO 6 M Q and ages of several millions up to 
several tens of million years, and sizes of parsecs, which yield 
support to an scenario where many instantaneous bursts oc- 
cur with a spread in ages. 



5 SN 2000FT AND OTHER VERY LUMINOUS 
RADIO SUPERNOVAE IN LIRGS 

SN 2000ft in NGC 7469 belongs, together with SN 2004ip 
in IRAS 18293-3413 (Mattila et al. 2007, Perez-Torres et al. 
2007) and SN 2008cs in IRAS 17138-1017 (Kankare et al. 
2008a, Perez- Torres et al. 2008, Kankare et al. 2008b), to the 
class of CCSNe that have been detected at both radio and 
optical/IR wavelengths in the nuclear regions of LIRGs. All 
of them have in common that are subject to a substantial ex- 
tinction, and show a high radio luminosity. A question arises: 
Is this simply a selection effect? The most luminous events 
at IR wavelengths are more easily detected by NIR searches 
that we are currently undertaking. The high IR luminos- 
ity can be expected to be powered by an interaction with 
a dense CSM, and the same scenario applies for high radio 
luminosity events. Thus, significantly less luminous events 
would simply remain undetected by current VLA observa- 
tions. Alternatively, it could be that such SNe are favoured 
in the nuclear regions of LIRGs, e.g. because of a top-heavy 
IMF, or a high metallicity, resulting in large mass loss rates 
for the SN progenitor stars, and therefore in a dense CSM. 
To get an unambiguous answer we need to keep searching 



for SNe in LIRGs at both IR and radio wavelengths so that 
we will have enough statistics in a few years from now. 



6 SUMMARY 

We have monitored the radio emission of the Luminous In- 
frared Galaxy (LIRG) NGC 7469 and its radio bright super- 
nova SN 2000ft from 1998 until 2006, using the VLA at 8.4 
GHz in A configuration. 

Our 8.4 GHz VLA monitoring of SN 2000ft for almost 
six years (t m 2100 days) has allowed us to show that its 
radio emission follows a rather steep decline (/3 ~ —2.02; 
Su oc t 13 ), which is typical of supernovae exploding in "nor- 
mal" environments. This result implies that the late time 
radio emission of SN 2000ft is still being powered by its in- 
teraction with the presupernova stellar wind, rather than 
by interaction with the interstellar medium (ISM). In fact, 
we find that the ram pressure of the presupernova wind is 
rather large p w w 7.6 x 10 -9 dyncm -2 , even at a super- 
nova age of t ^ 2127 days. This value is about three times 
larger than the inferred ISM pressure in the starburst re- 
gions of M 82 at a similar distance. Therefore, the pressure 
of the ISM where SN 2000ft exploded must be smaller than 
the ram pressure of the presupernova stellar wind, as oth- 
erwise we should have seen evidence of a flattening of the 
8.4 GHz flux density time decay. The swept-up mass by the 
supernova shock, after 2127 days of assumed free expansion 
in a steady, spherically symmetric wind (s — 2; p oc r" s ), is 
M sw w 0.29 Mq , and implies that the mass ejected in the ex- 
plosion of SN 2000ft must have been significantly larger. At 
this SN age, the shock has reached a distance r s h ~ 0.06 pc, 
and our 8.4 GHz VLA observations probe the interaction of 
the SN with material that was ejected by the progenitor star 
(via its stellar wind), about 5820 yr prior to its explosion. 
This wind material has a particle density n w ~ 5100 cm" 3 , 
which is similar to expected values of tiism in the dusty, 
dense environments of starbursts. 

We also searched in our radio images for recently ex- 
ploded core-collapse supernovae (CCSNe). Apart from SN 
2000ft (Lp C ak ~ 1.1 x 10 28 erg s" 1 Hz" 1 ), we found no 
evidence for RSNe more luminous than about (L poak w 
6.0 x 10 26 erg s _1 Hz" 1 ), suggesting that no other Type Iln 
SN has exploded since 2000 in the circumnuclear starburst 
of NGC 7469. 
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